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Summary. 1. We have investigated a physiological com- 
ponent of the gravitaxis of Paramecium using established 
mechanisms of ciliate mechanosensitivity. The horizon- 
tal, up and down swimming rates of cells, and the sedi- 
mentation of immobilized specimens were determined. 
Weak DC voltage gradients were applied to predeter- 
mine the Paramecium swimming direction. 

2. An observed steady swimming rate is the vector 
sum of active propulsion (P), a possible gravity-depen- 
dent change in swimming rate (A), and rate of sedimenta- 
tion (S). We approximated P from horizontal swimming. 
S was measured after cell immobilization. 

3. Theory predicts that the difference between the 
down and up swimming rates, divided by two, equals 
the sum of S and A. A is supposed to be the arithmetic 
mean of two subcomponents, J„ and Ap, from gravisti- 
mulation of the anterior and posterior cell ends, respec- 
tively. 

4. A negative value of A (0.038 irnn/s) was isolated 
with J„(0.070mm/s) subtracting from downward swim- 
ming, and Jp (0.005 mm/s) adding to upward propulsion. 
The data agree with one out of three possible ways of 
gravisensory transduction : outward deformation of the 
mechanically sensitive 'lower' soma membrane. We call 
the response a negative gravikinesis because both Ag and 
Ap antagonize sedimentation. 
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Abbreviations: A gravity-induced change in active propulsion; 
change in active propulsion as induced by anterior gravireceptor; 
Ap change in active propulsion as induced by posterior gravirecep- 
tor; D rate of downward swimming; D,- individual downward swim- 
ming rate; n number of data; P active propulsion of cell; $ mean 
swimming angle; p probability; R scalar value of mean swimming 
rate; Tq coefficient of orientation; r, coefficient of taxis response; 
S rate of sedimentation; ft individual angle of swimming direction; 
U rate of upward swimming; U, individual upward swimming rate; 
Vi individual swimming rate 
* To whom offprint requests should be sent 



Introduction 

Negative gravitaxis ^ in Paramecium is a spectacular up- 
ward orienting behaviour which has intrigued experi- 
menters and theorists since 100 years (Verworn 1889). 
The specific density of this cell is about 4% above that 
of common freshwater (Taneda 1987) which would in- 
duce slow sedimentation of cell populations in the ab- 
sence of compensating mechanisms. Oriented upward 
swimming of Paramecium was claimed to be effected by 
physical mechanisms (such as buoy-principle, hydrody- 
namics of small bodies, dislocation between centers of 
effort and gravity), or involving sensory transduction (re- 
viewed by Bean 1984; Haupt 1962; Machemer and De 
Peyer 1977). A large body of controversial and apparent- 
ly inconclusive data on gravitaxis in the early literature 
seems to have discouraged the search for guiding mecha- 
nisms for decades. Since cellular electrophysiology has 
become applicable to cihates, gravireception can be 
viewed in the light of mechanosensory transduction lo- 
cated in the soma membrane (Naitoh and Eckert 1969; 
see Machemer and Deitmer 1985). The bipolar organiza- 
tion of this mechanosensitivity favours a classical view 
of physiologically mediated gravitaxis, the statocyst hy- 
pothesis (Loeb 1897; Lyon 1905). Recent gravitaxis ex- 
periments in Paramecium and invertebrate larvae in fact 
suggest the existence of a gravity-induced cellular motor 
response (Mogami et al. 1988a, b; Baba et al. 1989). The 
present study takes a combined theoretical and experi- 
mental approach to cellular gravitaxis. Proceeding on 
established grounds of ciliary electromotor coupling, we 
postulate, that a gravity stimulus modulates the steady- 
state membrane potential and hence affects ciliary fre- 
quency, and the rate of forward swimming in Parame- 
cium. We have designed mass experiments and have tak- 
en extensive methodical steps to isolate and secure a 
gravity induced motor response. Our data prove the ex- 
istence of a negative gravikinesis in Paramecium and 



^ Traditionally, the term of 'geotaxis' was used in the literature. 
Because gravity is not restricted to terrestrial environments, 'gravi- 
taxis' has now been commonly adopted as a descriptive term 
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present a testable theoretical outline of transduction of 
the gravity stimulus. Preliminary data have been pub- 
lished in abstract form (Machemer 1989, 1990 a; Ma- 
chemer-Rohnisch 1989; Braucker and Machemer-Roh- 
nisch 1990). 

Theory 

Bipolar mechanosensitivity. The anterior and posterior 
ends of the cell soma of Paramecium are sites of high 
mechanosensitivity. Activation of anterior Ca mechano- 
receptor channels generates a depolarizing receptor po- 
tential, similar stimulation of posterior K mechanorecep- 
tor channels hyperpolarizes the membrane (Machemer 
and Ogura 1979). Simultaneous deformation of mechani- 
cally sensitive areas of the soma leads to summing of 
receptor potentials. Favourable cable conditions in Par- 
amecium spread the receptor potential involving virtually 
no losses, so that the motor response is activated simul- 
taneously in the cilia. The polarity and amplitude of 
the receptor potential determine the type of the ciliary 
motor response (Machemer 1988 a, b). Minor depolariza- 
tions decrease the rate of forward swimming, minor hy- 
perpolarizations slightly enhance forward swimming 
(Machemer and Sugjno 1989). 

We consider the possibility that the raised density 
of the cytosol (with respect to that of the medium) and 
its non-ideal viscous properties tend to deform the apical 
membrane inward (net force directed inward), and the 
'bottom' membrane outward (net force directed out- 
ward). Therefore, gravity may affect the anterior or pos- 
terior mechanoreceptors or both. We anticipate that ac- 
tivated anterior and/or posterior receptor conductances 
generate a tonic depolarizing, hyperpolarizing or no re- 
ceptor potential at all, depending on cell orientation and 
adequacy of stimulation. In principle, the receptor mem- 
brane may be responsive to one out of at least 3 types 
of gravity-induced mechanical load: (1) from both out- 
side and inside the cell, (2) from outside only, and (3) 
from inside only. 

Elements of gravity-induced swimming. Locomotion of a 
ciliated cell results from the vector sum of two forces: 
propulsion by ciliary activity (P) and sedimentation (5). 




Fig. lA-C. Determinants of the swimming rate of Paramecium in 
the gravity field (ant = anterior end of cell). Vectors of active pro- 
pulsion (P) and (downward) sedimentation (S) may be supple- 
mented by a gravity-induced kinetic component, A, which adds 
to, or subtracts from P. \. P+A may take any direction depending 
on cell orientation. R In upward swimming cells, the effective swim- 
ming rate is determined by direct subtraction of S from P + A ; 
during downward swimming, S adds to P + J. C. Sedimentation 
is not represented in a horizontal swimming trace; A is minimal 
or virtually absent (so that P is approximated) 



If the cell can sense the gravity stimulus and respond 
to it, a stimulus increment (or decrement) will add to 
propulsion (P+d) (Fig. lA). If the cell moves vertically, 
the contributions of S, A and P add numerically 
(Fig. 1 B). S does not affect horizontal displacement, and 
A approximates zero during horizontal locomotion; min- 
imal values of A may nevertheless persist under certain 
conditions (Fig. IC; see Discussion). 

Three alternative models of sensory-motor coupling under 
gravity 

Assuming that the cytosol acts as a heavy body of non- 
ideal fluid which tends to deform the top and bottom 
membrane, an inward pressure component exists across 
the upper membrane; the pressure is directed outward 

across the lower membrane. 

Model # 1 assumes that inward as well as outward defor- 
mation of the sensitive membrane is an adequate stimu- 
lus for gravitransduction, and that the electric membrane 
response is related to stimulus size (directed inward or 
outward) over the sensitive area of the membrane (Fig. 2 
left). During downward swimming (D) the speed of loco- 
motion is determined by propulsion (P), sedimentation 
(S), a gravity-induced hyperpolarizing conductance in- 
crease at the posterior cell end which induces an incre- 
ment in propulsion (zip), and a gravity-induced depolar- 
izing conductance increase at the anterior end, inducing 
a decrement in propulsion (/!„): 

Di = Pi4-zl4 + Zl„t + Si. (1) 
Similar conditions apply during upward swimming 
([/) (Fig. 2 middle): 

l/t=Pt+^pt+^ai + Si. (2) 

where and sedimentation act opposite to propulsion 
P and A p. Addition of Eqs. (1) and (2) give: 

Di-H/t = 2Si, or 

(D-U)/2 = S. 




Down Up Horizontal 



Fig. 2. Gravitransduction according to model # 1. Inward and out- 
ward deformation {thin arrows) of the mechanically sensitive soma 
membrane can elicit a gravireceptor response {+: depolarizing; 
— : hyperpolarizing). hft: During downward swimming (heavy ar- 
row), hyperpolarizing activation of posterior gravireceptors tends 
to raise the swimming rate (A pi), and depolarizing activation of 
anterior gravireceptors tends to decrease the swimming rate (AJ). 
The summed values of A^l and A„'\ add to propulsion and sedimen- 
tation (Eq. 1). Middle: Similar relationships apply to upward swim- 
ming (Eq. 2). Right: During horizontal swimming, the gravitational 
load on the anterior and posterior ends of the cell is minimal so 
that the net gravikinetic response approximates zero, and the active 
propulsion (P) is represented by the net swimming rate 
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Down Up Horizontal 



Fig. 3. Gravitransduction according to model # 2. Only inward de- 
formation of the sensitive anterior and posterior membrane can 
cause gravitransduction. In this model zlpj, and zl„J, share the posi- 
tive sign with sedimentation enhancing the downward swimming 
rates. Other definitions' corresponding to Fig. 2. See Eqs. (5) and 
(6) 
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Down Up Horizontal 



Fig. 4. Gravitransduction according to model #3. Only outward 
deformation of the sensitive anterior and posterior membrane can 
cause gravitransduction. In this model A^l and zl„t share the same 
sign, but unlike model # 2, they antagonize sedimentation during 
swimming. Other definitions corresponding to Fig, 2. See Eqs. (8) 
and (9) 



Model #1 predicts that the rate of sedimentation 
will be 50% of the difference between down and up swim- 
ming rates. Equation (3) is vahd also for the cases where 
the absolute amounts of Ap and A„ are different, or one 
subcomponent or both Ap and A^ are nil. 

During horizontal swimming the mechanical load on 
the sensitive anterior and posterior membranes will be 
at minimum (Fig. 2 right). In Paramecium, stimuli ap- 
plied to these areas polarize the membrane inversely so 
that we assume, as a first approximation, that Ap and 
Aa cancel out, and horizontal swimming approximates 
active propulsion (P). Then, subtraction of Eq. (2) from 
Eq.(l) leads to: 

iD + U)/2 = P. (4) 

Equation (4) states that horizontal swimming corre- 
sponds to the arithmetic mean of downward and upward 
swimming. Applicability of model # 1 to gravisensation 
in Paramecium can be tested by determinations of D, 
U, P and S, and using the Eqs. (3) and (4). 

Model #2 differs from model #1 in that it assumes 
gravisensory transduction only by inward deformation 
of the membrane (Fig. 3). Ap and A^ share the positive 
sign with sedimentation enhancing the downward swim- 
ming rates. Downward and upward swimming (Fig. 3, 
left, middle) are now defined as follows : 

Di=Pi + Api + Si, or Ap=D-P-S; (5) 
Ut = P'\ + A„i + Si, or A^ = P-U-S. (6) 

Addition of Eqs. (5) and (6), 

iD-U)/2 = S + {Ap+A„y2, (7) 

shows that with model # 2 the difference between D and 
U, divided by 2, exceeds sedimentation by the arithmetic 
mean of Ap and A„ which is equal to the value of A 
(Fig. 1). AppUcabihty of model #2 can be tested by mea- 
surements of D, U, P and S. Stimulus conditions during 
horizontal swimming correspond to those of model # 1 
so that the value of P is approximated. If the data show 
that {D-U)/2>S, the scalar values of Ap and A^ are 
derived from Eqs. (5) and (6). 



Model #5 assumes that gravisensory transduction acts 
by outward deformation of the sensory membrane 
(Fig. 4). Like in model #2, Ap and A„ have the same 
sign, but different from that model, they antagonize sedi- 
mentation during swimming. Downward and upward 
swimming are now defined as: 

Di=Pi + A„'[ + Si, or A„ = P-D + S; (8) 
U]=P'[ + Ap'[ + Si, or Ap=U-P + S. (9) 
Addition of Eqs. (8) and (9) give 

(D-Uy2^S-{Ap + AJ/2. (10) 

Equation (10) predicts that (D - {7)/2 < S.Hence, 
model # 3 differs from models # 1 and # 2, and can also 
be tested experimentally. Stimulus conditions during 
horizontal swimming are comparable to models # 1 and 
#2 allowing to approximate P. Additional determina- 
tions of D, U and S give the scalar values of Ap and 

Selection of the appropriate model. Statistically significant 
values of the half-difference between downward and up- 
ward swimming rates, and sedimentation, are the crucial 
experimental parameters determining which model is 
valid. If (D — U)/2 equals S, model # 1 is appropriate, 
or the swimming rate was unaffected by gravity. Mod- 
el # 1 does not resolve the fractional contributions of 
the anterior and posterior gravireceptors to motility. If 
(D — U)/2 is shown to be unequal to S, either model # 2 
(>S) or model #3 (<S) can apply. The sign and size 
of the arithmetic mean of Ap and A„ are generalized gra- 
vikinetic values, which may be expressed by A (Fig. 1): 

A=iAp + AJ/2. (11) 

Basic equation of vertical swimming. Paramecium cells can 
be induced to swim downward or upward by application 
of a vertical DC voltage gradient. The cells head toward 
the cathode irrespective of their location at the bottom 
or top of the fluid space (Fig. 5). Specificity of the mecha- 
noreceptor prevents interference of weak voltage stimuli. 
A low level of the applied voltage gradient also keeps 
potential-sensitive ciliary responses below saturation 
(Machemer 1988c). 
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Down 




Fig. 5. Induction of downward and upward swimming applying 
a vertical DC voltage gradient (0.6V/cm). During galvanotaxis, 
a cell swims forward heading toward the cathode. Simple switching 
of polarity of the voltage gradient reorients swimming from 'Down' 
to 'Up'. P, S and A add or subtract comparable to Fig. IB 



Figure 5 shows that the downward swimming rate 
is the vector sum of P, S and A (Eq. 12). The same applies 
to upward swimming with the direction of P reversed 
and the sign of A not reversed (Eq. 13; see models #2 
and #3). Summing of Eqs. (12) and (13) leads to Eq. (14), 
which expresses that the difference between downward 
(D) and upward swimming (C7), divided by 2, gives the 
sum of sedimentation and the gravity-induced mean kin- 
esis, A : 



Di=Pi + Sl + A (12) 

U'[ = P1 + Si + A (13) 
Di + J7T = 2(Si + zl), 

or {D-U)/2 = S + A. (14) 

Equation (14) represents Eqs. (3), (7) and (10) in a gener- 
alized form. The resulting sign and value of A identify 
the type of gravisensory transduction (models # 1 to #3) 
and determine the direction and size of gravikinesis. 



Materials and methods 

Culture. Wild-type of Paramecium caudatum were reared in Cero- 
phyl solution (Cerophyl Laboratories, Inc., Kansas City; 0.2% w/v 
cerophyl powder in aqua bidest., autoclaved), buffered at 7.0 pH 
by Sorensen buffer (1.8 mM Na2HPO4+0.2 mM KHjPOJ, bacter- 
ized with Aerobacter aerogenes, cultured at 22° C in 14 h/10 h light 
regimen, and harvested in early stationary phase after 2 days. 

Transfer to experimental solution. Cells were washed and collected 
after negative gravitactic acciunulation in a column of experimental 
solution (ImM CaClz + lmM KOH-|-1.6mM MOPS acid, 
pH 7.0). This washing procedure was repeated once, then the cells 
were infused into the experimental chamber (Fig. 6). Smooth 
transfer to the chamber of minimally agitated cells was achieved 
by 3 cm level difference of the solutions. 



A 'I 





Fig. 6A-C. Setup for recording graviresponses. A. 
Chamber (plexiglass; surface view) with dimensions in 
mm as indicated. The well includes a central part, 26 mm 
wide, for the experimental solution (entry: hi) and two 
peripheral sections, which include the electrodes (e; 
75 mm apart) and are filled with 1.5% agar in 
experimental solution through a bore {ha). A silicone 
rubber seal (s) surrounding the well is tightened by 8 
screws connecting cover and bottom of the chamber. B. 
Longitudinal section across chamber bottom and cover 
showing reduction to 1.6 mm height in the central part of 
the well. A 25-nm Oj and COj-permeable transparent foil 
intervenes between chamber bottom and cover (not 
shown). Omitted are also bores (3 mm in diameter) 
perforating the cover for access of air to the intercalated 
foil. C. Principle of recording vertical (continuous lines) 
and horizontal cell movement (dashed hues). A macro- 
objective (o) of the videocamera (c) focuses on the central 
area of the chamber (c), which receives oblique 
illumination (/) against a black background (f>) 
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Fig. 7. Variables affecting rates of horizontal swimming 
in the experimental chamber. The effects of culturing 
time (2 days: open symbols; 6 days: closed symbols), and 
time of ionic equilibration prior to incubation (up to 
14min: circles; 3h: triangles; 16h: squares) were tested 
to determine the appropriate time of cell incubation in 
the chamber before starting the experimental sequence. 
Swimming rates decreased during a period of 4 h and 
settled near 0.8 mm/s. The slow relaxation of swimming 
was largely independent of the time of incubation in the 
experimental solution suggesting that it was primarily 
determined by the O2— CO2 balance in the chamber. 
Motion affected by the closed chamber environment, 
was at steady-state after 4h of incubation 



Equilibration. A sample of about 200 cells was equilibrated and 
left unstirred for 4h in the chamber. During equilibration the 
chamber was kept in a water-saturated atmosphere at 22° C. An 
equilibration time of 4h is advised to secure full electrochemical 
equilibrium of the cell after transfer to a modified ionic environment 
(Oka et al. 1986; see Machemer 1990b). The same time was required 
for establishing a steady-state in — COj balance and hence swim- 



ming rate in the chamber environment (Fig. 7). Care was taken 
to fully isolate vibration from the chamber until the end of the 
experiment. 

Immobilization of cells. For determination of the sedimentation rate, 
samples of cells were exposed to a solution of 1 mM NiClj in experi- 
mental solution after the end of the active swimming period (Fig. 8). 
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Fig. 8. Experimental protocol for determination of the rates of 
swimming and sedimentation. The procedure is represented by time 
bars (read in rows from left to right). Recording starts with the 
chamber oriented horizontally (A) and then vertically (B, both unbi- 
ased swimming). The voltage gradient is turned on after 5 min to 
collect cells near the cathode (light arrow), followed by two 30 s 
sequences of swimming in opposite directions (C; heavy arrows). 
The chamber is then gradually reoriented by 90°, and the gradient 
is turned on to collect cells up; 3 pairs of 30 s down and 30 s 
up swimming are recorded (D). The chamber is reoriented horizon- 
tally, the gradient turned on to collect cells down, and a pair of 
horizontal swimming recorded (E). Further experiments recording 
unbiased horizontal (F) and vertical svsrimming (G). The chamber 



I 

is perfused for 5 min with 1 mM Ni^ * in experimental solution 
while the voltage gradient is on to save active cells from being 
swept out of the chamber. The cells are immobilized after further 
15 min of exposure to Ni^*. Sedimentation of the immobilized cells 
with the voltage gradient on is recorded during two time intervals 
of 4 min each (H). A second pair of 4-min sedimentations of the 
same preparation occurs with the gradient turned off (I). Seven 
parameters and circumstantial conditions were tested: 1. spatial 
symmetry of the chamber (A); 2. effects of time (A, F; B, G; C, E); 

3. effect of the voltage gradient on swimming rate (A, C; E, F); 

4. horizontal versus vertical swimming (C, D, E); 5. down versus 
up swimming (D); 6. completeness of <»11 immobilization (H, I); 
7. rate of sedimentation (H, I) 
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An amount of 8 ml of Ni^ "^-solution (corresponding to 5 chamber 
volumes) was slowly infused by gravity force from the lower side 
of the vertically oriented chamber. In order to reduce the loss of 
cells, a voltage gradient of 0.6 V/cm (cathode down) was applied 
for 5 min to accumulate active cells against the direction of infusion. 
At the end of the infusion period, 15 more min were allowed for 
completion of cell immobilization. According to microscopic tests, 
a minor fraction of cells (< 1%) shows residual ciliary activity after 
this time. Immobilized cells settle preferentially with their long axes 
vertical and their rear ends down (Kuznicki 1968; Roberts 1970), 
and this appUes to young cultures in the first line. In 
Ni^"^ -immobilized cells from 2-day cultures, Tancda et al. (1987) 
counted 66% of the specimens in vertical ( + 20°) sectors and 14% 
in horizontal ( + 30°) sectors. Roberts (1970) inferred a maximal 
velocity variation of 30% related to orientation for Paramecium. 
Because the majority of immobilized cells settles at vertical orienta- 
tion, the remaining fraction of specimens at oblique and horizontal 
orientations does not significantly shift the median of the sedimen- 
tation rate (see Fig. 13). This conclusion receives additional support 
from the low Reynolds number (< 10" ') of protozoan cells. 

Voltage gradient and galvanotaxis. A DC voltage gradient (0.6 V/cm) 
was applied during part of the swimming and sedimentation experi- 
ments. The gradient was stabilized using a constant current source 
supplying +150VDC for regulation. The current was calibrated 
so as to generate the predetermined voltage drop across the medi- 
um. Persistence of the current was monitored during the experi- 
ment. Generation of a voltage gradient by means of constant cur- 
rent overcomes time-dependent instabilities of the electrodes. The 
gradient of 0.6 V/cm (or 14 mV per cell length) corresponds to a 
maximal cathodal depolarization of 7mV (and maximal anodal 
hyperpolarization of the same size) for a cell oriented parallel to 
the field lines. Thus, the appUed gradient is a weak electric stimulus 
(Machemer 1988 b). 

Chamber. The experimental vessel (65 x 110 x 15 mm outer dimen- 
sion; Fig. 6) includes a fluid space of 41 x 85 x 1.6 mm or 5.6 cm^. 
This volume was aerated by diffusion using a gas-permeable 25 jim 
foil (Biofoil 25, Heraeus, Hanau) intercalated between the well and 
the perforated plexiglass cover (138 bores of 3 mm diameter). The 
aerated area of the fluid space corresponds to a 1.8 cm~^ surface-to- 
volume ratio equivalent to the upper surface of a 5.5 mm deep 
well in a circular Petri dish. Aeration strongly affects the swimming 
rate of the cells (Fig. 7). A 1.6 mm depth of the fluid space was 
chosen to minimize hydrodynamic wall effects (unstirred layer) on 
swimming rates and to secure focus for video recording of the 
cells. Chlorided silver wires were inserted at a distance of 75 mm 
from each other, in order to generate a linear voltage drop across 
the longer axis of the chamber. The electrodes were embedded 
in 1.5% agar in experimental solution. A central 8.5 x 11.5 mm field 
of the chamber was used for recording of swimming cells (for sedi- 
mentation: 4x5mm field). Prior to transfer of the cells, the 
chamber was filled with experimental solution to avoid the inclu- 
sion of air bubbles. The cells were laterally illuminated by infrared- 
filtered light. 

Recording. Cell movements were recorded using a CCD-videoca- 
mera (Panasonic FIO) with 25 Hz framing rate connected to a mac- 
ro-lens. Essential parameters such as date, time (10 ms resolution), 
horizontal or vertical orientation of chamber, and gradient polarity 
were superimposed on the video field for identification of the experi- 



Experimental schedule. Figure 8 explains the sequence of experimen- 
tal steps to record the swimming rates and sedimentation of a 
sample of about 100 cells. The schedule includes tests for: (1) ex- 
istence of a spatial asymmetry of the chamber (in terms of physical 
and/or chemical properties), which might affect the up and down 
swimming rate; (2) the effect of time on swimming rate; (3) the 
effect of the voltage gradient on swimming rate; (4) horizontal and 
vertical swimming rates; (5) down and up swimming rates; (6) com- 



pleteness of cell immobilization, and (7) rate of sedimentation. After 
turning on the voltage gradient of 0.6 V/cm, a 30 s period was 
allowed for cathodal accumulation of cells. The same gradient was 
applied (cathode up) during the sedimentation period to provide 
improved separation of remaining active cells from fully immobi- 
lized cells. Care was taken to avoid mechanical disturbances of 
the cells during changes between horizontal and vertical orientation 
of the chamber; a time of up to one minute was allowed for this 
90°-reorientation. 

Digital image analysis. Our processing of video records involves 
two automatic and one semi-automatic steps: 

1. Digitization of every 5th video-field (200 ms intervals; 512 x 256 
pixels; 256 intensity levels; RDDK-card, Bartscher, Eschwege); 
storage on hard disk together with time information. 

2. Apphcation to digitized fields of background-recognition algo- 
rithm, filtering (area-sensitive and common skeleton procedure) for 
isolation of cell shapes; superposition of 20 frames using colour 
sequence to show time. 

3. Marking of track start and end points by operator; automatic 
identification of cell boundaries, and monitoring of time sequence 
and steadiness of swimming speed by program; automatic calcula- 
tion of time intervals, length and angle of tracks, and swimming 
rates; automatic data storage. 

Statistics. For evaluation, we employ methods of non-parametric 
statistics because the distributions of swimming rates and angles 
of locomotion are unknown. Data samples are represented by medi- 
an; 95% confidence intervals describe distribution ranges. Experi- 
mental units (Fig. 8 A, F), (B, G), (C, E) were compared using the 
Mann- Whitney {/-test (two-tailed error probabihty of 5%; Sachs 
1984) to test for potential differences of the medians. Similar checks 
were done for the existence of a spatial bias of the experimental 
chamber (horizontal swimming rates in units C, F), and complete- 
ness of immobilization (units H, I). 

Model # 1 postulates the identity of the up and down swim- 
ming rates after compensation for sedimentation. Model #1 is 
tested according to Eq. (3) by comparison of the medians of individ- 
ual (Di — S) and (C7j+S) values. A statistically secured difference 
between these means rejects model # 1 and confirms the existence 
of A. 

Means of angles, and coefficients of orientation and taxis. The mean 
of directional data (individual angles /?,; « = number of data) is 
described by the mean angle * and the value of the mean vector 
R (Batschelet 1981). The coordinates of the 'centre of mass' (x, y) 
are calculated as: 

^ Z(cosfi) . ^. KsinjS,) . 




if x>0: arctan(y/x) 

if .x<0: 180° + arctan(j;/x) 

if jc = 0 : special further conditions. 

The orienting coefficient was defined as: 
ro = Rcos -P. 

Hence, ro=+l, if all cells are strictly oriented up (0°), ro= — 1, 
if all cells are oriented down (180°), and ?"o = 0, if all cells are oriented 
90° or 270°. 

For evaluation of the angular distribution of swimming rates, 
the frequency of swimming angles (p,) together with swimming rates 
(u,) must be taken into account. Thus, the values of R and $ were 
weighted by the swimming rate: 

^_ i:[(cosA)t>J . ^_ I[(sinA)fJ . 
(Definitions of R and * as shown above). 
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A taxis coefficient r, was defined ai 
-Rnc 



This coefficient combines tlie orientational response ( = mean angle 
of orientation) with the kinetic response of the cells. Hence, the 
coefficient represents the cellular taxis response. The value of r, 
equals —I, if all cells show negative (upward) gravitaxis; it is equal 
to +1, if all cells show positive (downward) gravitaxis, and the 
coefficient equals 0, if a gravitaxis response cannot be confirmed. 



Chamber symmetry test 

Physical and chemical properties of the fluid volume of 
the experimental chamber should not convey a bias to 
cellular swimming rates in the gravity field. The angular 
distribution of swimming directions and rates of the hori- 
zontally oriented chamber (see unit A in Fig. 8) were 
plotted in polar diagrams using 10° unit sectors (Fig. 9). 
Figure 9A shows that swimming between angles of 270° 
and 90° ('distal' direction corresponding to 'up' in a 
vertically oriented chamber) was equally frequent, and 
included similar angular distributions, as swimming be- 
tween 90° and 270° ('proximal', corresponding to 'down' 
in the vertical chamber). The mean angle (n = 899) was 
0°. The swimming velocities show an even distribution 
of directions (Fig. 9B). We express combined angular dis- 
tributions of frequency and swimming rates by a taxis 
coefficient ( + 0.09, horizontal chamber). This coefficient 
represents the negative y-value of the resulting mean ve- 
locity vector (Table 1). A coefficient value of -1-0.09 sug- 
gests that a taxis did not exist in horizontal swimming 
behaviour. 



Table 1. Coefficients related to gravitaxis. We calculate coefficients 
using the y-values (reversed sign) of firequency^weighted vectors 
from the data (Fig. 10, see Methods) or ratios of A over P (Table 2). 
The coefficient of gravitaxis incorporates both orientation and ki- 
netic components. Negative sign indicates effect antagonizing sedi- 
mentation. Contributions of orientation and kinesis are approxi- 
mated in fractions of the taxis coefficient. Coefficients of orientation 
and kinesis do not fully add to —0.23 because values include round- 
ing errors 



Experiment 


Coefficients of 




Orientation Kinesis 


Taxis 


Vertical, no V-grnd. 


-0.21(9/10) 


-0.23(1/1) 


Vertical, ^gradient 


-0.03'(1/10) 




Up, ^gradient 


-0.005''(2/100) 




Down, P^gradient 


-0.06 -=(3/10) 





° value of A divided by horizontal swimming rate (P) 
" value of Ap divided by P 
* value of A„ divided by P 



Time effects 

Several time-dependent variables may affect cellular lo- 
comotion: age of the Paramecium culture, length of the 
ionic equilibration period in experimental solution, and 
incubation time in the fluid volume of the chamber 
(Fig. 7). The culturing period (2 and 6 days since start 
of the culture) had no perceptible effect on the swimming 
rate (compare open symbols with closed ones), which 
stabilized near 0.8 mm/s after 4 h of incubation in the 
chamber. Ionic equilibration (prior to incubation) was 
not sufficient to achieve steady swimming rates; in any 
case it took 4 h of incubation in the chamber to stabilize 
swimming rates. The test (Fig. 7) suggests that a stable 




Fig. 9 A, B. Unbiased horizontal swimming (n = 899) as represented 
in polar distribution diagrams of frequency of swimming directions 
(A) and rates (B). The field was divided into 10° sectors. An angle 
of 0° corresponds to the orientation of 'distal' (equivalent to 'up' 
in the vertical swimming mode), and 180" correspond to 'proximal' 
(corresponding to 'down' in the vertical mode). A. Frequencies 
of orientations are the same in the proximal and distal halves of 



the diagram (resulting vector: 0.07 at an angle of 0°; no laterality!). 
Preponderance of laterally swimming cells reflects dimensions of 
the chamber fluid volume (Fig. 6) and screen dimensions for evalua- 
tions. B. Median swimming rates of sectors distribute evenly (result- 
ing vector: 0.1 mm/s, frequency-weighted angle of 230°). The taxis 
coefficient of -h0.09 suggests absence of directionality 
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physico-chemical environment, in particular a saturating 
O2 — CO2 gas exchange across the chamber membrane, 
exists after 4 h of incubation of cells in the chamber vol- 
ume (Fig. 6). 

The experimental protocol provides tests for unbi- 
ased swimming movement before and after recording of 
vertical swimming under the voltage gradient control 
(Fig. 8). No differences in the means of horizontal (A, F) 
and vertical swimming (B, G) were statistically secured, 
nor were the rates of horizontal swimming including a 
voltage gradient (C, E) affected by the intervening experi- 
ment (D). 



Effects of the voltage gradient on swimming rate 

The median rate of unbiased horizontal swimming was 
0.834 mm/s (n = 899). Applying the vohage gradient of 
0.6 V/cm raised the swimming rate by 36% (1.136 nmi/s, 
n = 1588; Table 2). Note that the velocity bias of the volt- 
age gradient cancels for determinations of the gravity- 
induced component A (Eq. 14). 



Vertical versus horizontal swimming 

In the absence of a voltage gradient, cells swam predomi- 
nantly upward confirming the existence of a negative 
gravitactic response (Fig. lOA). Angular velocity distri- 
butions ranged evenly between 0.8 and 0.9 mm/s 
(Fig. lOB; «=1038) with a median of 0.817 mm/s (Ta- 
ble 2). This median value resembles that of horizontal 
swimming (0.834 mm/s) excluding the possibility that 
gravity unilaterally raises or depresses the swimming 
rate. Frequency-related weighting of the resulting veloci- 
ty vector in Fig. 10 B (0.23 mm/s, 0°) might suggest that 
velocity tends to be reduced in downward swimming 
cells, and to rise during upward swimming in accordance 
with the negative gravitaxis coefficient ( — 0.23 ; gravitaxis 
coefficient in the horizontal state: +0.09). The reverse 
is true: cells heading upward between 270° and 90° were 
slower than cells swimming downward between 90° and 
270° (0.808 mm/s versus 1.045 mm/s), and this divergence 
grew, when the range of sectors was reduced to 10° 
(0.815 mm/s versus 1.045 mm/s; Table 2). The prelimi- 
nary analysis of horizontal and vertical swimming under- 



Table 2. Paramecium gravikinesis experimental results 



Movement 


^gradient 
applied? 




Displacement Rate [mm/s] 
median (confidence range') 


Horizontal 


no 


899 


0.834(0.813; 0.835) 


Vertical 




1038 


0.817(0.799; 0.835) 


Up (270° -90°) 




674 of 1038 


0.808(0.784; 0.829) 


Down (90° -270°) 




364 of 1038 


0.843(0.810; 0.881) 


Up (355°- 5°) 




49 of 1038 


0.815(0.725; 0.937) 


Down (175°- 185°) 




19 of 1038 


1.045(0.635; 1.246) 


Up 


yes 


1552 


1.057(1.034; 1.079) 


Down 


yes 


1512 


1.150(1.132; 1.180) 


Horizontal 


yes 


1588 


1.136(1.120; 1.154) 


Sedimentation 


no/yes 


1362 


0.084(0.081; 0.086) 


A 






0.038 








0.070 








0.005 



' 95% confidence interval of the median value 




Fig. 10 A, B. Vertical swimming in 
the absence of a voltage gradient 
(n=1038). A. Angular distribution 
of frequency shows a preference 
of upward swimming. The 
resulting vector (0.21, 0°) 
corresponds to upward 
orientation of swimming. 
B. The angular swimming rates 
are evenly distributed between 0.8 
and 0.9 mm/s (a few cells headed 
downward near 1 mm/s). The 
resulting vector (0.23 mm/s, 
frequency-weighted angle: 0°) 
shows a moderate upward bias of 
swimming velocities (taxis 
coefficient: -0.23) precisely 
inverse to the gravity vector 
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scores the need for improved separation of the orienta- 
tional and kinetic components during gravitactic swim- 
ming of Paramecium. 



Down and up swimming rates 

In the presence of a voltage gradient, values of vertical 
and horizontal swimming rates (horizontal: 1.136 mm/s; 
up: 1.057 mm/s; down: 1.150 mm/s; Table 2) confirm the 
shift in velocity obtained without voltage gradient: 
downward swimming rates were slightly larger, upward 
swimming rates were smaller than those during horizon- 
tal swimming. Polarograms of frequency and velocity 
illustrate that the voltage gradient acted to restrict swim- 
ming directions (Fig. IIA, C) and rates (Fig. UB, D) to 
approximately +15° from the vertical axis. Swimming 
rates outside the +15° sector were rare). 

With n>1500, the Up and Down swimming rates 
were clearly separated from each other (P< 0.001%; 
Fig. 12 A; Table 2); the overall range of velocities corre- 
sponds to that seen in horizontal swimming (Fig. 12B). 
Because sedimentation subtracts from upward swim- 
ming, and adds to downward swimming, the observed 



difference in net up and down rates might be effected 
by the gravity vector. We have tested this possibility 
adding the sedimentation rate (0.084 mm/s; Table 2) to 
the Up-value, and subtracting the same rate from the 
Down-value. The resulting values of active propulsion 
(Up: 1.141 mm/s; Down: 1.066 mm/s) are still statisti- 
cally significant (P< 0.001%), and are thus a strong indi- 
cation for a physiologically mediated cellular response 
to gravity. 



Completeness of immobilization, and absence 
of electrophoresis 

Cells, which are exposed to 1 mM NiClj for 20 min, be- 
come immobilized with their morphological properties 
unchanged. Preliminary tests have shown that some cells 
retain a fraction of active cilia, which might convey an 
undesirable bias to the sedimentation data. We have ap- 
plied a voltage gradient of 0.6 V/cm to the immobilized 
cells (cathode up) expecting that the remaining minor 
population of active cells tends to swim against sedimen- 
tation thereby distorting the sedimentation distribution 
histogram (Fig. 13). Sedimentation rates of the same cells 




Fig. IIA-D. Vertical swimming 
(A, B: upward, n=1552; 
C, D: downward, n= 1512) in the 
presence of a voltage gradient 
(0.6 V/cm; cathode and swimming 
direction coincide). A. Upward 
swimming was largely limited to 
±15° of the 0°-line (resulting 
vector: 0.95, 0°). B. Angular 
distribution of swimming rates 
with the cathode up. The 
resulting vector was 1.02 mm/s at 
a frequency-weighted angle of 0°. 
Note that lateral and even 
downward swimming did occur in 
a very small number of cells. The 
taxis coefficient of —0.95 
indicates a near perfect upward 
orientation, which is based, in 
this case, on galvanotaxis. 
C. Frequency pattern of 
downward swimming angles is 
inverse to that of upward 
swimming (resulting vector: 0.87, 
170°). D. Distribution of 
swimming rates with the cathode 
down is inverse to (B). Resulting 
vector: 0.85 mm/s, 180°; taxis 
coefficient: -1-0.85 
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FREQUENCY 




SWIMMING RATE [mm/s] 



DOWN (n = 1512) CZ] UP (n = 1552) 



FREQUENCY 




SWIMMING RATE [mm/s] 



B H PROXIMAL (n = 729) DISTAL ( n = 859) 

Fig. 12A, B. Histogram representing swimming rates of equilibrat- 
ed Paramecium applying a 0.6V/cm voltage gradient. A. Down 
(«= 1512) and Up rates {n= 1552) cover a wide range of individual 
velocities. Note that a positive shift of the Down-rates with respect 
to the Up-rates applies to all velocities. B. Horizontal swimming 
rates. Proximal (n = 729) corresponds to Down, and Distal (n = 859) 
to Up. Comparison with (A) shows that swimming rate distribu- 
tions are unrelated to gravity 



H. Machemer at al. : Gravikinesis in Paramecium 

in the presence and absence of the voltage gradient were 
identical, however, suggesting that cell immobilization 
was sufficiently complete for safe assessment of sedimen- 
tation. Immobilized cells were not moved in voltage gra- 
dients up to 2V/cm. This excludes, together with the 
observations mentioned above, the existence of a signifi- 
cant amount of electrophoresis in immobiUzed and Hve 
paramecia during their exposure to the weak voltage 
gradient. 

Rate of sedimentation. Cells sedimented in a wide range 
of velocities 0.035 to 0.345 mm/s (Fig. 13). The median 
of the distribution was 0.084 mm/s (n=1362, Table 2). 
Statistical analysis shows the median of sedimentation 
rate (of our 2-day cell populations in experimental solu- 
tion) is stable down to n = 100; even with «=25, the 
median varied between 0.067 and 0.097 mm/s only. 

Isolation of the gravity-induced kinetic component. With 
the values of D (downward rate), U (upward rate) and 
S (sedimentation rate) determined, Eq. (14) yields an ab- 
solute value of A of 0.038 mm/s which eliminates mod- 
el # 1 in addition to the statistical test mentioned pre- 
viously. The sign of A is negative excluding model #2 
and confirming model #3. A negative value of 
0.038 mm/s means that the mean gravikinetic component 
antagonizes sedimentation. 

Model #3 assumes, by electrophysiological reason- 
ing, that the value of A is the arithmetic mean of two 
independent subcomponents, A„ and Ap (equation 10). 
Application of equations (8) and (9), that is, inclusion 
of the value of P (unbiased horizontal propulsion) into 
calculations, isolates the scalar values of A„ 
= 0.070 mm/s, and of Jp= 0.005 mm/s (Table 2). These 
data suggest that downward swimming cells slow down 
by so as to largely compensate the effect of sedimenta- 
tion (0.084 mm/s), and upward swimming cells slightly 
speed up by Ap. The value of Ap is small and may be 
non-existent. 
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Discussion 

The present paper is the first quantitative and statisti- 
cally secured account of a physiological component in 
gravitaxis of Paramecium. The existence of a gravikinetic 
component was suggested by previous experiments (Mo- 
gami et al. 1988 a; Baba et al. 1989), and was postulated 
on theoretical grounds (Machemer 1989). Mogami et al. 
(1988a) observed a minor reduction of the swimming 
rate of Paramecium during free fall (n = 9), and this reduc- 
tion was more pronounced in cells swimming upward 
(n = 6). A tendency in Paramecium for active compensa- 
tion of sedimentation was also inferred from the behav- 
iour of slowly swimming cells at 1 G, and of fast swim- 
ming cells in the presence of hypergravity ( > 3 G ; Baba 
et al. 1989). Anoxic conditions raised the swimming ve- 
locity and improved gravitaxis in Paramecium (Hem- 
mersbach-Krause et al. 1990); it remains to be deter- 
mined, if gravitaxis profited from the kinetic response 
or was specifically affected by anoxia. A strategy for iso- 
lation of a gravikinetic component according to Eq. (14) 
was proposed by Machemer (1989), and preliminary re- 
sults confirming the existence of this component were 
presented by Machemer-Rohnisch (1989). 

We attempt to further secure the results of the present 
paper by combination of the mass-cell approach with 
an effort to control a variety of circumstantial parame- 
ters such as culture conditions, temperature, symmetry 
of the chamber, ionic equilibration, gas exchange, me- 
chanical disturbances, electrophoresis. Our experimental 
design, subtraction of upward and downward swimming 
rates, further eliminates potential biasing factors such 
as effects of the electric DC field on swimming, and bio- 
convection. We sometimes observed phenomena resem- 
bling convection in sedimenting aggregates of Ni-treated 
immobile cells. Trajectories of these records were dis- 
carded. Convection was not observed in isolated immo- 
bile or active paramecia. If bioconvection nevertheless 
did occur in cathodally swimming cells, this artefact is 
cancelled due to the subtraction procedure. 

Our data analysis suggests that the motor responses 
of individual paramecia to gravity may vary. This is no 
surprise since in large populations of unsolicited cells, 
horizontal as well as vertical swimming rates cover a 
wide range (Fig. 12), and so do individual sedimentation 
rates (Fig. 13). According to experimental experience, de- 
polarizing and hyperpolarizing mechanoreceptor cur- 
rents scatter in individual specimens (Ogura and Ma- 
chemer 1980; Machemer-Rohnisch and Machemer 
1984); similar conditions may apply to gravireception. 
The presumed individual variability in electric mem- 
brane steady-state conditions (size of leakage conduc- 
tance) and sensitivity (type and number of channels 
available for gravireception) can explain, why some cells 
swam actively downward (Fig. lOA) exceeding rates of 
upward swimmers by up to 40% (Fig. lOB). Individual 
variability of experimental cells cautions the experi- 
menter to draw generalizing conclusions from too small 
populations of cells. Our present data are based on ap- 
proximately 1000 different cells (some of our readings 
(«) may include the same cell over different intervals). 
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Statistical analysis of these data shows that 400 cells 
would lead to reliable conclusions; for sedimentation 100 
specimens would suffice. 

We were not prepared to see the major effect of gra- 
vikinesis referred to depolarizing activation (of the pre- 
sumed anterior gravireceptor) during downward swim- 
ming, whereas hyperpolarizing activation (of the pre- 
sumed posterior gravireceptor) was minimal. The rele- 
vance of this difference may depend on the goodness 
of the value of unbiased active propulsion (P; Eqs. 8 
and 9). We have assumed that P is reasonably approxi- 
mated during horizontal swimming in the 1 — G environ- 
ment. The real value of P needs to be determined in 
the virtual absence of the gravity vector (such as during 
free fall or in orbit). If P was overestimated by only 
3%, both A„ and Ap approximate A; the latter value 
remains unaffected by P (Eq. 14). 

The sizes of A^ and Ap might have been altered by 
the voltage gradient which generated, for a 230 \im cell, 
a (maximal) 7 mV hyperpolarization at the posterior cell 
apex during upward swimming, and a (maximal) 7 mV 
depolarization at the anterior cell end during downward 
swimming. A voltage dependence of mechanoreceptor 
conductances (gjcgca) has been shown in Stylonychia: 
the K-dependent hyperpolarizing conductance (g^) de- 
creased at potentials more positive than the resting po- 
tential (F„r), and the Ca-dependent depolarizing conduc- 
tance (gca) decreased with potentials negative to V^, (see 
Machemer and Deitmer 1985). Extrapolation of these 
data to a hypothetical minor F-dependence of receptor 
conductances in Paramecium (about 1% change per mV 
shift in potential) would not apply to our voltage gra- 
dient of 0.6 V/cm ; this size of gradient shifts the sensitive 
membrane to potentials which do not affect the conduc- 
tance of the mechanoreceptor. A weak desensitizing ef- 
fect of pre-hyperpolarization on the posterior mechano- 
receptor (and, by inference, on Ap) cannot be excluded 
a priori because hyperpolarizing conditioning of the 
membrane slightly reduces the K driving force and the 
size of the hyperpolarizing receptor potential. V„r ap- 
proximates — 30mV, and — 90mV for our experi- 
mental solution (Ogura and Machemer 1980). 

Using the calculated sizes of J„ and /d^ at face value, 
they make sense : starved cells which are little responsive 
to an upward reorienting mechanism (Taneda et al. 
1987), or cells having unusual distributions of Ca and 
Mg phosphorous crystals (Kaneshiro and Rope 1989) 
and/or lipids (Cole et al. 1990), may obviate fast sedimen- 
tation by means of reduction of their downward swim- 
ming rate. Inversely, the significance of the kinetic re- 
sponse, A„, decreases with a rising proportion of orienta- 
tion in gravitaxis suggesting that these mechanisms can 
partly substitute each other. 

Our results show that gravitaxis in Paramecium is 
a complex behaviour including orientational as well as 
kinetic responses. We employ coefficients for quantita- 
tive description of gravitaxis (Table 1) showing that this 
response is not all or none. Our taxis coefficient of — 0.23 
suggests, however, that a cell population has a slow net 
upward drift. The proportion of the orientational com- 
ponent in establishing gravitaxis, as compared to the 
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kinetic component, appears to be in the order of 9/10 
versus 1/10. The kinetic component potentially rises to 
a proportion of about 1/4 of gravitaxis in downward 
swimming cells. 

In establishing the existence of a physiologically me- 
diated response to graviy, we are not claiming that gravi- 
taxis is based on gravisensation only. We would rather 
assume that the continuing controversy on either physi- 
cal or physiological foundations of gravitaxis will even- 
tually lead ad absurdum because both views are correct: 
(1) postulating gravisensation by interaction of cyto- 
plasm and membrane (statocyst hypothesis; Loeb 1897; 
Lyon 1905; Koehler 1922), and (2) explaining gravitaxis 
in Paramecium on mechanical principles (Verworn 1889; 
Roberts 1970; Winet and Jahn 1974; Fukui and Asai 
1985). 
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